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Abstract. We discuss the production of charmonium in nuclear collisions within the 
framework of the statistical hadronization model. We demonstrate that the model 
reproduces very well the availble data at RHIC. We provide predictions for the LHC 
Qjj ■ energy where, dependently on the charm production cross section, a dramatically 

different behaviour of charmonium production as a function of centrality might be 
expected. We discuss also the case in elementary collisions, where clearly the statistical 
model does not reproduce the measurements. 



1. Introduction 

Charmonium production is considered, since the original proposal more than 20 years 
ago about its suppression in a Quark-Gluon Plasma (QGP) pQ, an important probe 
to determine the degree of deconfinement reached in the fireball produced in ultra- 
relativistic nucleus-nucleus collisions. In the original scenario of J/ip suppression via 
Debye screening p] it is assumed that the charmonia are rapidly formed in initial hard 
collisions but are subsequently destroyed in the QGP (see recent summaries in [21 E]). 

In recent publications jl] we have demonstrated that the data on J/-0 and ip' 
production in nucleus-nucleus collisions at the SPS (^snn ~ 17 GeV) and RHIC 
(v / sw r /v=200 GeV) energies can be well described within the statistical hadronization 
model proposed in [6] and have provided predictions for the LHC energy ( A /svA r =5.5 
TeV) and for energies close to threshold (s/snn ~ 6 GeV). Because the mass of heavy 
quarks exceeds the transition temperature of the QCD phase transition by a factor of 
almost 8, heavy flavor hadron production cannot be described in a purely thermal 
approach. It was, however, realized in [6] that charmonium and charmed hadron 
production can be well described by assuming that all charm quarks are produced 
in initial, hard collisions while charmed hadron and charmonium production takes 
place exclusively at the phase boundary with statistical weights calculated in a thermal 
approach. For a recent review of this statistical hadronization approach see [5]. An 
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important element is thermal equilibration, at least near the critical temperature, T c , 
which we believe can be achieved efficiently for charm only in the QGP. 

An important question is whether the statistical behavior is a unique feature of 
high energy nucleus-nucleus collisions or whether it is also encountered in elementary 
collisions, where finite size effects will obscure a possible phase transition. In early 
analyses (see ref. [7]), it is indeed argued that hadron production in e + e~ and pp is 
thermal in nature. Furthermore, such analyses of hadron multiplicities (for recent results 
see [HI El [TO] and refs. therein) yield also temperature values in the range of 160-170 
MeV. 

To establish the uniqueness of the J/i[) probe for the diagnosis of a QGP in nucleus- 
nucleus collisions it is important to understand whether similar thermal features as 
observed in nucleus-nucleus collisions are also at work in heavy-flavor hadron production 
in elementary collisions. We have recently demonstrated [TTJ that there are crucial 
differences in elementary collisions compared to nucleus-nucleus. 

2. The statistical model in elementary collisions 

For the study of hadron production in e + e~ collisions we employ the canonical statistical 
model described in [HI [TTJ (see also [9j H2]). For the present study, we perform 
calculations for two cases: i) a 2-jet initial state which carries the quantum numbers of 
the 5 flavors, with the relative abundance of the five flavors in one jet and corresponding 
antiflavor in the other jet taken form the measurements at the Z° resonance quoted in 
[T3] . These relative abundances (17.6% for uu and cc and 21.6% for dd, ss and bb) 
are thus external input values, unrelated with the thermal model, ii) a purely thermal 
ansatz , i.e. a 2-jet initial state characterized by vanishing quantum numbers in each 
jet. 

In Fig. [JJwe show a comparison of data [13] and model prediction for charmed and 
bottom hadron yields in e + e~ annihilations at a/s=91 GeV. For the model we have used 
the parameter set: T=170 MeV, V^=16 fm 3 and 7 S =0.66, which represents the best fit 
of multiplicities of hadrons with lighter quarks [H] . 

We first note that the calculation employing the 5-flavor scheme is in very good 
agreement with the data, as demonstrated by the good \ 2 P er degree of freedom between 
the model and the data (excluding the T and prompt J/V) of 21.7/16 (34/18 when 
including all species) (see also ref. [9]). Despite this overall agreement, the exceptions are 
significant: the T meson yield is underpredicted by the model by 17 orders of magnitude, 
while the prompt J/i/j yield [TT] is underpredicted by almost 2 orders of magnitude. 
Obviously, the production of quarkonia is expected to be strongly suppressed in the 
statistical model. The disagreement is a consequence of the separate hadronization of 
the c and c quarks. The measured prompt J/ip production in Z° decays (into hadrons) 
is about 3xl0~ 4 [H]. The thermal model predicts a prompt yield for J/ip of 4.1xl0 -6 
(1.6xl0 -7 for ip' and 4.3xl0 -7 for Xd), identically for the two calculation schemes. 
Whenever the model seems to describe the yields of charmonia the measured yields are 
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Figure 1. Multiplicities of 
hadrons with charm and bottom 
quarks in e + e~ collisions com- 
pared to the thermal model cal- 
culations for two cases: i) the 5- 
flavor jet scheme (thick lines) and 
ii) no (net) flavor jet scheme (thin 
lines with diamonds for the charm 
sector and triangles for bottom). 
Note, for case ii) the large factors 
used to scale the model calcula- 
tions to fit in the plotting range. 
The data are from the compilation 
published by the Particle Data 
Group (PDG) [13 . The prompt 
J/tp measurement J/^ pr is from 
the L3 experiment [14] . 

dominated by the feed down from bottom hadrons and the agreement only reflects the 
agreement seen for the open bottom hadrons and their branching ratios to charmonia, 
properly considered in the model. 

The calculation employing a purely thermal ansatz underpredicts all the 
measurements by many oders of magnitude, while for the light quark sector the 
differences between calculations with a pure thermal model and with the 5-flavor quark- 
antiquark scheme were found to be small [8]. The strangeness suppression factor, which 
for the present results only enters in the calculation of the yields of D s and B s mesons, 
appears to have no counterpart in the heavy quark sector. This reflects the fact that 
a negligible number of c and b quarks are formed in the fragmentation process. In this 
case, the thermal weights describe the distribution of the initial quarks into heavy flavor 
hadrons. Chemical equilibration is not required in this process. 

For the case of hadron production in elementary hadronic collisions we employ 
the canonical realization of the thermal model (151 HI [16]. For the description of the 
relative production cross sections of heavy flavored hadrons, the energy dependence of 
the temperature parameter is the only model input, which is taken in a parametrized 
form from the fits of (ii,<i,s)-hadron abundancies in central nucleus-nucleus collisions 
|17j . For cm. energies beyond 10 GeV per nucleon pair in nucleus-nucleus collisions a 
limiting temperature T; im =164±5 MeV is reached. Recent fits of hadron yields in pp 
collisions [10] give very similar values, independent of anergy. 

In Fig. [2] we show the model comparison to data for the relative production cross 
section of ip' and J/ip charmonialjl The measurements in pA and pp(p) collisions are 
above the model values by about a factor 4 (corresponding to 10 experimental standard 

| The charm production cross section, which is an important model input parameter for the calculations 
of absolute yields, cancels out for this ratio. 
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Figure 2. Production cross section 
of if)' relative to J/ if). The data for 
pA collisions are from the compilation 
by Maltoni et al. [19 ; the points 
for elementary collisions are from the 
PHENIX experiment at RHIC (20] and 
from the CDF experiment at Tevatron 
[22] (see text); the data point for Pb+Pb 
collisions at the SPS energy is from the 
NA50 experiment [18]. The average 
value of the pA and pp(p) measurements 
with the corresponding error (see text) is 
represented by the shaded box. The band 
denotes statistical model calculations for 
the temperature parametrization from 
heavy-ion fits 17] (T ;im =164 MeV) with 
±5 MeV error. 

deviations; the average value of the measurements is 0.137±0.009, with a \ 2 P er degree 
of freedom of 0.88). The relative production cross sections of charmonium states, as 
are observed in all measurements in hadronic collisions cannot be described in the 
thermal approach. The temperature needed to explain the data would be 300 MeV, 
well above the Hagedorn limiting temperature, which is about 200 MeV. This is in 
sharp contrast to the (only currently existing) measurement in central nucleus-nucleus 
collisions, performed at the SPS by the NA50 experiment [IB] , which is well described. 
We note that the pA data exhibit a constant ip'/J/ip production ratio as a function 
of energy. In the model, the value is determined only by the temperature and this is 
reflected in the slight decrease of the ratio towards low energies. A constant value, also 
up to the LHC energies, is predicted beyond ^snn —20 GeV. 

The measurements reported in Fig. [^demonstrate that the relative production cross 
section ifj'/J/i/j is identical in pA and in pp(p) collisions, implying no visible influence 
of the cold nuclear medium. Note that the ratio for the Tevatron energy was derived 
from the CDF measurements of J/ip [21] and ip' [22] and is for transverse momentum 
p t >1.25 GeV/c (we have extrapolated the ip' measurement from 2 GeV/c down to 1.25 
GeV/c). 

3. The statistical hadronization model in nucleus-nucleus collisions 

The model has the following input parameters: i) characteristics at chemical freeze-out: 
temperature, T, baryochemical potential, //&, and volume corresponding to one unit of 
rapidity V& y= i, extracted from thermal fits of non-charmed hadrons [IT]; ii) the charm 
production cross section in pp collisions, taken either from NLO pQCD calculations 
[23] 124"] or from experiment [25], used to calculate the number of directly produced 
cc pairs N^ r . The heavy quark (charm or bottom) balance equation [6] (including 
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is used to determine the fugacity factor g c (J„ are modified Bessel functions; iV*^ and 
A^g are the numbers of open and hidden charm hadrons, respectively, in the volume 
Vaj/=i, computed from their grand-canonical densities). 

A comprehensive set of model predictions as a function of energy, from the charm 
production threshold up to the LHC, is presented in Fig. [3] for central Au-Au collisions 
(N part =350). The left panel shows our predictions for the energy dependence of 
midrapidity yields for various open and hidden charm hadrons. The most striking 
behavior is observed for the production of A+ baryons: their yield relative to other 
charmed hadrons rises significantly towards lower energies. In our approach this is 
caused by the increase in baryochemical potential towards lower energies (coupled with 
the charm neutrality condition). A similar behavior is seen for the baryon. The 
relative production yields of D-mesons depend on their quark content and depend on 
energy only around threshold. These results emphasize the importance of measuring, 
at low energies, in addition to D-mesons, also the yield of charmed baryons to get a 
complete measure of the total charm production cross section. 
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Figure 3. Left panel: energy dependence of the production yield at midrapidity for 
open and hidden charm hadrons (with the energy dependence of da cS /dy as in [4]). 
Right panel: ratios calculated from midrapidity yields. 



One of the motivations for the study of charm production at low energies 
was the expectation [26| [27] to provide, by a measurement of D-meson production 
near threshold, information on their possible in-medium modification near the phase 
boundary. However, the cross section a c5 is governed by the mass of the charm 
quark m c ~ 1.3 GeV, which is much larger than any soft QCD scale such as Aq C d- 
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Therefore we expect no medium effects on this quantity. The much later formed D- 
mesons, or other charmed hadrons, may well change their mass in the hot medium. 
Whatever the medium effects may be, they can, because of the charm conservation, 
&cc = k( a D + oa c + ...) + (& Vc + <Jj/ip + ■■■), expressed in Eq. ([1]), in first order only lead to 
a redistribution of charm quarks [I]. In contrast, as a consequence of in- medium masses 
of open charm hadrons, the yields of charmonia can vary by up to 40-50% and this is 
more prominent at threshold energies [I]. 

In the right panel of Fig. [3] we show the energy dependence of ratios of mid-rapidity 
yields of open and hidden charm hadrons and pions. The relative production of charmed 
hadrons to pions exhibits a monotonic dependence on energy, unlike equivalent ratios 
involving strange hadrons [IT]. The different bahavior compared to the strangeness 
sector is mainly due to the strong canonical suppression of open charm at low energies 
(determined in turn by strong decrease of the charm production cross section at low 
energies). A non- monotonic behavior is observed in the production ratio of J/ip to D°, 
determined by the interplay of the strong canonical suppression of open charm at low 
energies and the gradual onset as a function of energy of the quadratic term in Eq. (Tj[|). 
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Figure 4. Rapidity dependence of R J Jt f° r two centrality classes. The data from the 
PHENIX experiment [25J (symbols with errors) are compared to calculations (lines, 
see text). The shaded area corresponds to calculations for the lower limit of the charm 
cross section as measured by PHENIX [25] . with our shadowing scenario. 

In Fig. H]we present the rapidity dependence of the nuclear modification factor R^a 
(we use experimental data for the pp reference). While earlier pi] we have compared 
data [28] to our model predictions for the pQCD charm production cross section [23] 
(dashed line in Fig. H]), we use here as alternative the charm cross section as measured 
by PHENIX in pp collisions [25] and consider in addition shadowing for Au-Au collsions 
as extracted from recent dAu data [29] ; we assume that the deviation of R J Ja u from 
unity is entirely due to shadowing. Also in this case, our model describes the observed 
suppression and its rapidity dependence. The maximum of Raa & t midrapidity is in our 
model due to the enhanced generation of charmonium around mid-rapidity, determined 
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by the rapidity dependence of the charm production cross section. In this sense, the 
above result constitutes strong evidence for the statistical generation of J/ip at chemical 
freeze-out. 




Npart Npart 

Figure 5. Centrality dependence of R J /^ for RHIC and LHC energies (left panel) 
and of the J/tp rapidity density at LHC relative to the number of initially produced 
cc pairs (right panel, curves labelled by the da C c/dy) at midrapidity. 

The centrality dependence of i?^f at midrapidity is shown in the left panel of 
Fig. |5j Our calculations approach the value in pp collisions around N part =50, which 
corresponds to an assumed minimal volume for the creation of QGP of 400 fm 3 [4] . The 
model reproduces very well the decreasing trend versus centrality seen in the RHIC data 
[28]- Note that, in our model, the centrality dependence of the nuclear modification 
factor arises entirely as a consequence of the still rather moderate rapidity density 
of initially produced charm quark pairs at RHIC (dA^c/dy=1.6 for central collisions, 
using the FONLL charm production cross section [23]) leading to a marked centrality 
dependence of the canonical suppression factor for open charm hadrons. At the much 
higher LHC energy the charm production cross section (including shadowing in PbPb 
collisions [21]) is expected to be about an order of magnitude larger. As a result, the 
opposite trend as a function of centrality is predicted, with R^t exceeding unity for 
central collisions. A significantly larger enhancement of about a factor of 2 is obtained if 
the charm production cross section is two times larger than presently assumed, as seen 
in the right panel of Fig. El where we show the J / ip production relative to the number 
of initially produced cc pairs. 

4. Conclusions 

We have shown that the statistical hadronization model describes well the measured 
decrease with centrality and the rapidity dependence of R J Ja a ^ RHIC energy. 



Heavy quark(onium) at LHC: the statistical hadronization case 



8 



Importantly, well described in central heavy-ion collisions at SPS energy, where we 
have the only measurement to date, is the ratio ip'/J/ip. In contrast, in elementary, 
e + e~ and in pp and p(7r)-nucleus, collisions the charmonium measurements cannot be 
explained within the statistical model. This underlines the relevance of the statistical 
features in nucleus-nucleus collisions. Extrapolation to LHC energy leads, contrary to 
the observations at RHIC, to R^a increasing with collision centrality. If observed, such 
a behavior will provide a dramatic proof of charmonium as the ultimate observable to 
delineate the phase boundary of QCD matter. 

K.R. acknowledges support from the Alexander von Humboldt Foundation. 
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